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I. INTRODUCTION 

It is believed that in two-dimensions (2D) even small disorder localizes all electron states.0 
First hint for the absence of the mobility threshold in 2D came from the calculation of 
the weak localization correction, 5a(u), to the conductivity at finite frequency to. It was 
showni that in the limit k F l ^> 1 one has 5a(u)/ao = ^ In |co>r|, where k F is the Fermi 
momentum, I is the mean free path, and r is the elastic scattering time; o"o = (e 2 /h)kpl 
is the Drude conductivity. The logarithmic singularity in 5a (u) indicates that the zero- 
temperature (T = 0) conductivity is metallic only if u ^> ujq, where 

co ~ -e-***' (1) 

T 

is the characteristic frequency which marks the crossover to the exponentially localized 
regime.S From the finite size, L, correction to the zero-frequency conductivity, which is of 
the order of \n(L/l), one can estimate the localization length £ as§ 

e = iexpf|Mj. (2) 

Therefore, starting from metal, one cannot describe the T = behavior of the conductivity 
at frequencies lj ~ u . The adequate language for this region would be the language of 
the localized states. Within this language the conductivity originates from the transitions 
between the localized states induced by an external a.c. field. 

Let us briefly remind the corresponding derivation of cr(u), carried out by MottH in the 
strongly localized regime. The Hamiltonian for strongly localized electrons reads 

H = Y^ e i c l c i + %(4 c j + C ) C i)i ( 3 ) 

i ij 

where £j is the energy of a localized state centered at r = rj and Jjj is the overlap integral 
which falls off exponentially with distance: 1^ = Joexp(— |r^ — Tj\/a), where a is the size of 
the wave-function of the localized electron. The general expression for the conductivity is 
given by the Kubo formula 
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*M = ^EIM>r g ^. n (4) 

where (i|#|j) is the matrix element of x calculated from the exact eigenstates |i) and |j) of 
the Hamiltonian @) with energies E- x and Ej, huj\j = E\ — E- } and n- x is the occupation number 
of the state |i). 

The dissipative conductivity, Recr(u;), is determined by the pairs of states with = uj. 
At huj <C Iq the spatial separation between bare states \i) and \j) is much larger than a. 
Then in the calculation of eigenstates of a resonant pair one should take into account the 
overlap iy within this pair only and neglect the overlap with all the other localized states. 
This gives 

|i> = ^|i>+^|j>, |j> = ^K>+^p|j>- (5) 
The corresponging energies are 

= (q + ej)/2 ± T/2, T = [fa - ej) 2 + = hu iy (6) 

Using (|5|) the matrix element in (H) takes the form 

= (xi - x^Iij/T. (7) 

The contribution to Kea(uj) from pairs with shoulder r can be presented as 

Rea(^, r) = ^^J*^ F(hu, r), (8) 

where F(hu>, r) is the density of pairs with shoulder r and excitation energy tku [here J(r) = 
Ioe~ r ^ a ). The density F(hu, r) is determined by the condition that the pair is singly occupied 
(with energies on the opposite sides from the Fermi level). Then we have 

F{h», r) = g> J de x J de 2 6 + *±*) - *±^) 5 [ yfe - e 2 )* + AP(r) - hu 

ggW (9) 
A /(^) 2 -4/ 2 (r)' 

where (7 is the density of states and 9(x) is the step-function. Substituting (Q) into (|8]) and 
integrating over r we get 



Rea(u) = e -{2ix 2 g 2 huo) / , U =. (10) 



For /icj <C Jo the main contribution to the integral comes from r ~ r u = ahx(2Io/tiu) ^> a 
and one obtains 

ReaH = V2~Tr 2 j(g 2 ah 2 uj 2 rl). (11) 



Although we cannot show it explicitly, we argue below that the Mott expression ( |iT|) 
for Kea(uj) is valid also for the Anderson insulator with kpl 3> 1. When applying ( |TT|) 
to the Anderson insulator one should replace a by £ from (0) and substitute g = m/2irh 2 
(m is the electron mass). However, the question remains: what is the magnitude of Jo in 
this case? A plausible estimate can be obtained in the spirit of the Thouless picture of 
localization.!^ By definition, I represents the splitting of energy levels of two neighboring 
localized states (with centers at distance ~ £). The estimate for I Q emerges if one equates 
this splitting to the mean energy spacing for localized states centered within the area ~ £ 2 , 
so that Iq ~ l/g£ 2 (see also Ref. |5]). 

Note that in ID case a similar argument leads to Iq ~ l/pi^, gi being the ID density 
of states. Important is that in ID the Kubo formula can be evaluated exactly§ resulting in 
the ID version of the Mott formula, from which one can recover the above estimate for Jo- 
Such a mapping was first established by Shklovskii and Efros.0 

Using Eq. ( PD one can formally evaluate Re<r(u;) for hu > 2J . In this case the main 
contribution to the integral comes from r ~ a and we get 

Rea(u) = ^(J Va 4 ). (12) 

Certainly, the presentation of J(r) in the form Joe _r//a makes sense only for r a This 



means that the numerical coefficient in Eq. (|12|) is not reliable. 



Clearly, at large frequencies u ^> u the conductivity of the Anderson insulator should 
have the Drude form. The fact that Eq. ([12]) , calculated for strongly localized electrons, is 
also frequency independent allows us to assume that the description of a.c. transport in the 



Anderson insulator based on Hamiltonian (^) is accurate within a numerical coefficient. In 
other words, we assume that despite the complex structure of the electron wave-functions 
in the Anderson insulator the energy dependence of the matrix elements calculated between 
these functions is still given by (|7|). 

Eq. flT2| ) provides yet another way to estimate Iq in the limit kpl ^> 1. Namely, it matches 
do if we take Jo ~ {kpl) 1 ^ 2 / g^ 2 ■ We see that the dependence of Iq on £ in both estimates is 
the same; the extra factor (kpl) 1 ^ 2 presumably can be accounted for the \n(kpl) corrections 
to the exponent of £. 

There is also another argument in favor of the above estimate for Iq. The frequency 
dependence of Rea(u) in the Anderson insulator becomes strong for u <C Uq, whereas in the 
picture of strongly localized electrons the demarcation frequency is u> ~ Iq/%. Equating I 
to Tiujq we get Iq = kpl/g!; 2 , with another extra factor kpl. 

The simplified description of the Anderson insulator based on the Hamiltonian ([3D allows 
one to include into consideration the Coulomb correlations (i.e., the correlations in the 
occupation numbers of the localized states caused by electron-electron interactions) using 
the ideas first spelled out in Refs. §],[7||| This is the main goal of the present paper. We 
study the effect of Coulomb correlations on both Recx(u;) and Imcr(w). The most drastic 
conclusion we come to is that due to modification of lma((j) by the Coulomb correlations a 
system of localized electrons can support surface plasmons within a certain frequency range. 
We also show that these plasmons cause an additional structure in the behavior of Rea(u) 

at UJ > OJq ~ Iq/%. 

The paper is organized as follows. In Section |TJ we analyze the polarizability of the 
localized system in the absence of Coulomb correlations. In Section |TJ we introduce the 
Coulomb correlations and find the dispersion law for the surface plasmons. In Section |IT1| 
we study the corrections to the dispersion law due to the resonant scattering of plasmons 
by pairs of localized states. In Section [TV] we calculate the plasmon contribution to Recr(co>). 
Section [V] concludes the paper. 
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II. POL ARIZ ABILITY IN THE ABSENCE OF COULOMB CORRELATIONS 



In this section we demonstrate that without Coulomb correlations the 2D Anderson 
insulator cannot support surface plasmon. 

In the framework of the linear response theory the dispersion law of a plasmon, u(q), is 
determined from the condition^ 

l=v(q)ReP(uj,q) J (13) 

where v(q) is the Fourier component of the electron-electron interaction v(r) and P(u,q) 
is the polarization operator. Within a standard approach P(u,q) is calculated for non- 
interacting electrons described by Hamiltonian (^): 

P K5 ) = i ?Ki |^i j)P _-^_. (14) 

In the absence of disorder the eigenstates |i) and |j) are the plain waves so that the matrix 
element in ( JH| ) reduces to the delta-function 8(i— j — q). Then evaluating P(u,q) and 
substituting it into ( |T3"| ) together with 2D Coulomb interaction v(q) = 2-ne 2 / nq (k is the 
dielectric constant) yields the surface plasmon with the well-known dispersion law 



u{q) 



27rne 2 q 



rriK 



1/2 



(15) 



where n is the 2D concentration of electrons. The plasmon mode is undamped if q < uj/vp 
where vp = (AimY^Ti/m is the Fermi velocity (for larger q the Landau damping leads to a 
finite IhlP). The latter condition can be rewritten as q < 1/2clb where clb = h 2 K/me 2 is the 
Bohr radius. 

In the case of a strong disorder the eigenstates i and j in ( |14D are the localized states. 
The polarization operator (|14|) can be evaluated in a way similar to that employed in the 
Introduction for calculation of a.c. conductivity. For small q the matrix element in (|14D can 
be evaluated using (0): 

(i|e^ r |j) = zqr/(r)/r. (16) 



Then the contribution to P(w, q) from the pairs with the shoulder r [density V(uj, q, r) of 
the polarization operator] can be presented as 

^.'■^WMjCf^ (17) 

where T and F(T,r) is given by Eqs. (^) and (|S[), respectively (the integral is understood as 
principal part). Substituting F(T,r) into (|17|) and integrating over T we obtain 

r 



poo 

ReV(u, q, r) = 2# 2 (qr) 2 / 2 (r) / dT 

J2I(r) 



2/(0 ^r 2 - 4/2(0 [(Ro;) 2 -r 2 ] 
7T(7 2 (qr) 2 / 2 (r) 



for hu < 21 (r) 



= 0, for huj > 2/(r). (18) 

We see that ReP(ti;, q, r) is either negative or zero, so that Eq. ( |13D cannot be satisfied. 

In fact, the result ([Tj|) is almost obvious. Indeed, the imaginary part of polarization 
operator density, ImP(w, q, r), at small q differs from Rea(u, r) in (§) by a factor e 2 u/q 2 , so 
it follows from (§) and @ that 

1mP(u,q,r) = ^£M2!ilL, for hu>2I(r), 

y/{hu) 2 - 4/ 2 (r) 

= 0, for ^ < 2/(r). (19) 

Since ReP(uj, q, r) and lmV(u>, q, r) are connected via the Kramers-Kronig relation, the form 
(|19|) immediately follows from (p~8|) . 



III. COULOMB CORRELATIONS AND SURFACE PLASMON 

In the previous section the polarization operator was evaluated using the pair density 
(§) which was derived for non-interacting electrons. As it was first pointed out by Efrosi, 
interactions modify strongly the density of singly-occupied pairs. The underlying physics 
is the following. A pair can be singly-occupied even if both energy states reside below the 
Fermi level. The right condition for the pair to be singly-occupied is that the addition of 
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a second electron (which interacts with the first one) is energetically unfavorable. Such a 
Coulomb correlations effectively enhance the density of "soft" pairs (i.e., the pairs with small 
excitation energy T). 

Our goal is to apply the latter argument, which was presented for strongly localized 
system, to the Anderson insulator with large £. In order to do so we will adopt two assump- 
tions: 

(i) The interactions do not change the localization radius £. 

(ii) The estimate for the overlap integral, Io ~ l/g^ 2 , is unchanged in the presence of 
interactions. 

In other words, we assume that switching on the interactions leads to the Coulomb shifts 
of the eigenenergies but does not affect the wave functions. Note that assumptions (z) and 



(ii) contradict those made in Refs. |TT| and 0, respectively. 

As we will see below, the relevant pairs would be those with shoulder r ~ £. In other 
words, the relevant transitions shift the position of electron by ~ £. To establish the form of 
the density of singly occupied pairs F(T,r) with such a shoulder one can argue as follows. 
An isolated region of a size £ can be viewed as a small metallic granule. The transfer of 
an additional electron into this granule leads to the charging energy U = e 2 /2C, where C 
is the capacitance of a granule. In other words, the levels in the granule get shifted^ by 
an amount ~ U. Consider now two neighboring granules and assume that U ^> Iq. Due 
to aforementioned charging effect the highest occupied levels in the two granules typically 
differ by ~ U. Let for concreteness the highest occupied level in the first granule be higher 
by U than in the second one. Then the sought singly occupied pair with frequency V can 
be composed from the top occupied states in the first granule (these states should belong to 
the energy interval T + U measured from the highest occupied level) and the empty states 
in the second granule. Then the density of pairs F(T, r) (with r ~ £) can be estimated as 
g 2 (T + U). With the energy splitting taken into account it can be written as 

/r 2 -4/ 2 (r) 



Then at U = we return to Certainly, our consideration, based on artificial arrang- 
ing the localized states into the granules, provides only the order of magnitude estimate 
of F(T,r). In particular, the numerical coefficient in ( pop cannot be found from such a 
consideration. Our choice of numerical coefficient in (^) provides matching with a similar 
expression for strongly localized regime .010 

With the pair density (|20|) we can now easily evaluate ReP(uj,q,r). Calculating the 
integral over T in (|17|) yields 



ReV(u, g, r) 



ng 2 (qr) 2 I 2 (r) 



/4/ 2 (r) - (hu) 2 
U 2g 2 (qr) 2 P(r) 

hu 



2 U 
71 hu 



arctan 



hu 



for hu < 21 (r) 



: In 



(hu) 



AP(r) 



(hu) 2 -AP(r) + hu 



2I(r) 



4/ 2 (r) - (hu) 2 . 

for hu > 2I(r). (21) 



The expression fl2T|) for ReV(u, q, r) can be also obtained, using the Kramers-Kronig rela- 
tions, from lmV(u, q, r) which has a simple form 



lmV(u, g, r) 



7r(qr) 2 / 2 (r) , 



(huf 



7rg 2 (qr) 2 I 2 (r) ( Z7 \ 
■ , KH 1 y = 1 + — , for hu > 21 (r) 

^(M 2 - 4/ 2 (r) V W 



0, for hu < 2I(r). 



(22) 



Note that as in flT9|), lmV(u, q, r) ^ only for fro; > 2I(r). 

We see that the enhancement in the density of pairs with small V leads to a positive sign 
of ReP(c<j, g, r) for frc<j > 21 (r). This is our main observation. 

In order to obtain ReP(u, q) one should integrate ReV(u, q, r) over r. For hu > 2/ the 



main contribution to this integral comes from r ~ £ [like in derivation of Eq. (12)] and we 
get 



_ p( , 3vr g 2 UP g 2 e . fftwN 

ReP( "' g) = T^M^ ln W 



(23) 



For generality we will assume that there is also a gate at a distance d from the plane of 2D 
electrons. In this case the Fourier component of electron-electron interaction has the form 



v(q) 



2ix& 
nq 



1 - e~ 2(ld 



(24) 



If the gate is close to the electron plane, that is qd <C 1, we can expand the exponent in fl24|) 
and get v(q) = 47re 2 d/K. If d <C £ then the capacitance C reduces to the capacitance of two 
disks with area S = £ 2 separated by a distance d, so that C = nS/Aitd and consequently 
U = 4:Tre 2 d/ k,£, 2 . With these v(q) and U, after substituting ( p3|) into the plasmon equation 
(|TB|), we obtain 

'-MS^Kt)- (25) 

Certainly, the numerical coefficient in (|25|) should not be taken seriously. According to the 
assumtion (ii), I ~ g£ 2 . Then Eq. (^) yields the following dispersion law for the surface 
plasmon 

■m-HSMB- (26) 

We see that the dispersion law is close to acoustic. 

Let us establish the frequency range for the surface plasmon with dispersion law (f26|). 
The validity of expansion (|TED, g£ <C 1, implies that U ^> TiujIyT 1 / 2 {Tiuj / Iq) . On the other 
hand, %u > 2J . Then the frequency range for plasmon is 2I Q < %u < U . The nesessary 
condition for this range to be wide is U Iq. The ratio U/Iq — 47re 2 c//£ 2 / with J Q ~ l/g£, 2 
can be presented as 8ir 2 d/aB, where Ob = Ti 2 n/me 2 is the Bohr raduis. Thus, the condition 
d ^> cib insures that the plasmon equation (|13D has a solution within a wide frequency range. 
If d < a,B, the screening of Coulomb interaction by the gate is strong and the number of soft 
pairs is not sufficient to change the sign of ReP(u;, q). 

A similar condition can be obtained from the analysis of ImP(ci;, q), which can be derived 
by integration of Eq. ( |22|) over r 

Im ^') = XW ( ^ + c/) ~(^ ( ^ + c,) - (27) 

The origin of ImP(w, q) is the interaction of a plasmon with "resonant" pairs having exci- 
tation energy uo. In fact, lmP(u, q) describes the resonant scattering of a plasmon by a pair 
of localized states. One can introduce a mean free path, I, associated with such a scattering 
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and obtain ql ~ ln 1 ^ 2 (/?,co'//o)/(l + hu/U). Thus, the condition 2I < hou < U reduces to the 
condition ql > 1. 

In the absense of gate [gc/ ^> 1 in Eq. ([24]) ] one should take {7 = e 2 /ft£, so that [//Jo ~ 
£/o,b ^ 1- Then after a simple algebra we obtain 

■m-HSMS- (28) 

It can be shown that the above analysis of the validity applies in this case as well and leads 
to the same frequency range 2 Jo ^ Tiuj ^ U. Within this range we again have g£ <C 1. 

Thus, one should use Eq. (|26|) for qd < 1 and Eq. (^|) for qd > 1. On the other hand, 
the magnitude of U depends on the ratio d/^. Note that for d > £ one still can have qd < 1. 
In this case the dispersion law is given by Eq. (|26|) with U = e 2 



IV. RENORMALIZATION OF THE PLASMON DISPERSION LAW 

In the previous section, when calculating the polarization operator, we took into account 
the Coulomb correlations within a pair, but neglected the effect of polarization of surround- 



ing pairs on a given pair. On the other hand, by averaging of the polarization operator 
over frequencies of pairs T and their shoulders r we have effectively replaced the localized 
system by a medium. The average polarization of this medium gave rise to a plasmon mode. 
Within this procedure the "feedback" from surrounding pairs reduces to the interaction of 
a given pair with plasmons. In the present section we study the renormalization of the 
plasmon spectrum due to this effect. 

Generally, the plasmon excitation is defined as a pole in the density-density correlation 
function, H(u>, q, q'), which is related to the polarization operator P(uj, q, q') by the Dyson 
equation 

n(a;, q, q') = P(u, q, q') + J 7^pP(^ qiM<Zi)II(u;, q 1; q'). (29) 
Before averaging, both P(a>, q, q') and IT(u;,q, q') depend on two momentum variables q 
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and q'. The approximation we made above reduces to replacing of P(o>, q, q') by its average 
P(o>, q), so that the solution of fl29|) takes the form 



n(w,g) = P }\$ v (30) 

1 l-v(q)P(u,q) K } 



Then the pole of H(uj, q) is determined by the plasmon equation flT3|). 



As a next step, we took for P(uj,q) its expression (|14D for non-interacting electrons, 
which represents a sum of polarizations of pairs, 

p( w ,,) = iEK^U>l'fl.M. W = STpt^M). (31) 

U J 

and performed the summation neglecting correlations between the pairs but with the 
Coulomb correlations within a pair included. 

The renormalized Py(&>) for a given pair can be obtained from the following procedure. 
The function H(u, q) has a diagrammatic presentation in a form of a series of bubbles (ij), 
corresponding to Py, connected by the Coulomb interaction lines (see Fig. la). First, we 
arrange into a single block the sum over all combinations of bubbles which appear between 
two bubbles (ij) (see Fig. lb). Then we replace this block by its average, so that the result 
can be presented as two bubbles (ij) connected by a plasmon propagator Tl(u>, q) (see Fig. lb) 
(there is also an extra factor v(q) in each vertex). The renormalized bubble (ij) can be then 
obtained by summing up the series, consisting from the bubbles (ij), connected by plasmon 
lines (see Fig. lc). The resulting expression for Py(o>) reads 

(32) 

with 

3uM = / ^l(i|e iqr |j)|V(g)n(^,g). (33) 
Finally, replacing Py(o?) in ( |3"ID by Py(o>) we obtain 

P^rt^^'^' ^ + ^-^-^M - (34) 
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The equations (p0[), (|33|), and (|34]) form a closed system which determines II(a;,q) and, 
correspondingly, the renormalized dispersion law of a plasmon in a self-consistent way The 
approximation made in order to get the closed system [replacement of a block by a sought 
function II (a>, q)] is known as the effective-medium approximation. 

The analysis of the system (^,^,^) reveals that the renormalization of the plasmon 
dispersion law is weak. Namely, the appearence of the term i?y(o;) in the denominator 
of ( |3~4|) has the physical meaning that a pair (ij) acquires a finite life-time ry due to the 
interaction with plasmons. We will show that this life-time is long, i.e., 1/ry <C cuy. It 
can be readily seen that the difference between the renormalized polarization ReP(co>, q) and 
ReP(a;, q) originates from resonant pairs with (u;— c^y) ~ 1/ry and r ~ £ (note that the main 
term, ReP(u, q), is determined by the entire interval cjy ~ uj). If we neglect the dependence 
of the matrix element in fl34|) on ojy then the renormalization correction to KeP(uj, q) would 
be identically zero. A finite correction results from a slight asymmetry of the matrix element 
within the narrow interval (uj — cuy) ~ V r y- Then the relative magnitude of the correction 
is of the order of 1/^Ty with Ty calculated for a pair with = uj. Expecting Ty to be long, 
we can calculate it by substituting the nonrenormalized dispersion law of a plasmon into 
Eq. fl3"3p. Performing the integration we obtain 

R^u) = - ~ huj ^-J ^-J In wrth gate, 

~hw^jj-\ (^JJ^J ^ n_3 ^"/ — ^' without gate. (35) 

Since both ratios, hu/U and Iq/U are small, the correction to the dispersion law, 
8q(uj)/q(uj) ~ 1/uJTjj, is negligible. 

In the next section we will see that the corresponding renormalization of Recr(u;) is much 
larger than the renormalization of the dispersion law. 
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V. RENORMALIZATION OF THE REAL PART OF THE CONDUCTIVITY 

As we have seen in the previous section, the renormalization of the polarization operator 



where the summation is performed over the singly occupied pairs ij and in this way the 
Coulomb correlations are taken into account. In the limit ry — > oo and for u > 2Iq one 
should use the pair density F(T,r), given by Eq. (|20|) , in order to perform the summation. 
The result is determined by resonant pairs with u = cjy = T/h: 



We see that within the frequency interval 2J < Tluj < U the real part of the conductivity 
exceeds the Drude value due to Coulomb correlations. 

When calculating the correction to Recr(u;) caused by the finite value of Ty it is important 
to realize that ft/ry is maximal for soft pairs with small u;y. This is because the matrix 
element (i|e iqr |j) is proportional to l/u;y. In the previous Section this was not important 
since the correction to ReP(u), q) came from the resonant pairs only. Here, however, we have 

Im(u + LU-jj + V r ij) _1 V r ij ^ l/^U' so so ^ P a i rs §i ye the main contribution to 

the correction 5Kea(uj). Assuming cuy <C to we can present this correction in the form 



results in an appearance of ih/r-^ in the denominator of Eq. ([34]). Correspondingly, the 
renormalized expression (|j) for ~Rea(u) can now be rewritten as 





(37) 



Recr(a;) 




E 



(^ij) 2 Tij 




(38) 



where we have substituted h/r^ = i?y from Eq. (|33|) . The sum over pairs is again evaluated 
with the pair density ( p0|) 




(39) 
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The main contribution to the integral over T comes from the lower limit T ~ Then the 
integral over r is again determined by r ~ £, so that the sum (^) appears to be ~ I^g 2 !] . 
As in the previous section, the value of the integral in ( p9|) takes different values in the 
presence and in the absence of a gate. Finally we obtain 

~ (if) (x) ln 3 (^)' wi thout gate. (40) 

Comparing ( [35]) to (|40|) we see that both corrections are of the same order (I /U) 2 at hoo ~ J . 
However the correction to Rea is much bigger at ^> Jo. This reveals a new mechanism 
of absorption of a.c. field: by resonant excitation of plasmons. More precisely, the field 
polarizes the soft pairs (with cuy ~ /o/^) and the induced polarization excites the plasmon 
waves. Therefore, the energy of a.c. field is effectively absorbed by plasmons. The rapid 
increase of 5a with u is caused by the number (phase volume) of plasmons which absorb 
the field. 

Note that with correction 5a(u) the total conductivity Recr(u;) exhibits a rather com- 
plicated behavior. For hu < 2Iq the conductivity increases with 10. Then it passes 
through a maximum at uo ~ lo/h and falls off with uj according to Eq. (0). How- 
ever at Tiuj ~ £/(/ /[/) 1//4 (with gate) and hto ~ U(Iq/U) 1 I & (without gate) we have 
<5Recr(co>) ~ Recr(u;) and the conductivity starts rising again. On the other hand, the expres- 
sion for 5Kea was derived assuming that it is small. Therefore in the region 5Kea > Rea(cj) 
the renormalization of <j(oj) by plasmons is strong. In this case one cannot calculate ry using 



the bare polarization operator. The full analysis of the system (|30| , |33| , |34]) in this frequency 
range is out of the scope of the present paper. 



VI. CONCLUSION 



In the present paper we argue that the wave of electric field can propagate along the 
surface of the 2D Anderson insulator. The field originates from the density fluctuations of 
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localized electrons. One should distinguish this wave from the usual plasmon in an ideal 
2D gas with the dispersion law given by Eq. fll5|) : (i) the derivation of (|T5| ) implies that 
uj ^> 1/r while we predict the existence of a plasmon at much lower frequencies uj ^ lo/^- 
The minimal frequencies 1/r and /o/^ differ by a factor exp(7r/ci?Z); (ii) the plasmon ( |T5|) 
results from the solution of Eq. ( |T3"D with polarization operator calculated for free electrons. 
Within this approximation Eq. (|l~3"D has no solutions for localized electrons. The solution 
appears only if one takes into account the Coulomb correlations in the occupation numbers 
of the localized states. 

The obvious consequence of the existence of a plasmon excitation is that localized elec- 
trons "feel" a fluctuating electric field S(u,q) with the spectral density 

n f( x,2\ jf_ q 2 v 2 (q)ImP(uj,q) 

) 27Te 2 [l-v{q)ReP{uJ,q)} 2 +[v{q)lmP(uJ,q)] 2, 1 ' 

where P(u;,g) is the polarization operator fl34|). The spectral density (f4~T|) has a peak at 

uj = uj(q) corresponding to the dispersion law of a plasmon, which is different with and 

without gate. 

The basic assumption of our theory is that in the presence of interactions the Anderson 
insulator is characterized by two energy scales: Iq - the spacing between energy levels in 
the area £ 2 and the charging energy U which is the Coulomb interaction of two localized 
electrons separated by a distance ~ £. This energy modifies the pair density F(uj,r). Such 
a modification results in an enhancementi^B of the dissipative conductivity Reo"(a;). If we 
denote as p(u) the matrix element of r calculated for a pair with frequency u then 

Recr(cu) oc u p 2 (uj) F (u , r w ) . (42) 

As it was discussed in the Introduction, it is plausible to assume that for the Anderson 
insulator in the absence of interactions p{uj) behaves as p{oj) ~ r u for uj < I /h and p{uj) ~ 
^Iq/uj for uj ^ Io/h. Since the pair density in the absence of interactions is proportional to 
uj we reproduce Eqs. ([□]) and (|12"D: 

Rea(co') oc uj 2 , for uj < Iq/Ti, (43a) 
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Recr(co>) = const, for uj ^ Io/h, (43b) 



Where we neglected a logarithmic factor in ( |43aj ) . We have assumed that the frequency 



dependence of p(uj) remains the same in the presence of interactions.^ At the same time, 
F(u;, £) is changed drastically by interactions: F(uj, £) = const foruJ<U /h and F(u>,£,) oc uj 
for uj > U/h. As a result, we get the following frequency dependence of Rea(u) in the 
presence of interactions [see also 



Rea(u) oc ou, for u < Io/h, (44a) 

Rea(u) oc for I /h < uj < U/h, (44b) 

Recr(co') = const, for uj > U/h, (44c) 

This simple analysis shows that Rea(uj) exhibits a maximum at uj ~ /o/^- ( see a ls° the end 
of the previous section). In fact, this maximum is intimately related to the existence of a 
plasmon. Indeed, for Eq. (|T3| ) to have solutions we need a positive sign of ReP(c<j, q), which 
is equivalent to a positive lma{uj). At the same time, lmo~(uj) can be obtained from Rea(uj) 
using the Kramers- Kronig relation. In order to trace how the modification of Re<r(u;) by 
Coulomb correlations leads to the change of sign of lma(uj), we can interpolate the frequency 
dependence of p(uj) as p{ui) oc \{hui) 2 + 4/q] -1 / 2 and the frequency dependence of F(u,£) as 
F(uj, £) oc {huj + U). Then we have 



^WtM' ^ (45) 



where G = -^-^ is the conductance. It can be easily seen that Eq. (|45|) reproduces correctly 
all the limiting cases (0) and ([44]) both with and without charging effect. The imaginary 



part of the conductance calculated from (|45|) has a simple form 



lmG(w) = (L)^ + 4 - u > °' (46) 

In Fig. 2 we have plotted R,eG(uj) and ImG(o;) for different ratios U/2I Q . We see that in the 
absence of Coulomb correlations (U = 0), lmG(uj) is strictly negative (in our calculation in 
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Section [TJ it turns to zero for u > la/%). However, at finite U the change of sign occurs at 
to = (2I /h) exp(iiI /U). For U ^> I Q this frequency is just 2I /h. 



As a final remark, let us outline the difference between our approach and that of Ref. 16 



In Ref. [16] the authors addressed electron-electron interactions in the strongly localized 
regime when the localization radius is much smaller than the interpair separation. In this 
case singly-occupied pairs can be considered as point-like dipoles. The dipole moment 
induced by an external field £oe~ tuJt can be written as 

rf^r^+E^i, (47) 

Ij^k 

where £ k is a component v of the electric field caused by polarization of a dipole I which 
acts on the dipole k (summation over repeating indices v is implied). The polarizability a% v 
of a dipole k has the form 

ay _ 2e 2 PkPk^k _ £ 2 a v p (ao\ 

ak ~lT^^I-n PkPk ^ P ^ (48) 

ij 

where p k and u k are correspondingly the matrix element and the frequency of the dipole k. 
Here Py is given by Eq. (|3l|) with the states i and j making up the dipole k. To keep the 
discussion simple we will assume that the polarizability is isotropic, i.e., a^ u = The 
field £u can, in turn, be expressed through p fc as 



£, 



(49) 

R=Rt R; 



dRv R 3 

where R& is the position of the dipole k. 

Upon substituting (^) into ([17]) we obtain an infinite system of linear equations. It can 
be easily shown that iterating this system leads to the renormalization of polarizabilities of 
dipoles a k 

d R v \( d R»\ 

H . (50b) 

R=Rfc— R; 



ak 



8R» R 3 / 1 dR v R 3 
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Here is the self-energy. Fleishman and AndersonS anylized this self-energy using the 
arguments similar to those put forward by Anderson0 when he demonstrated the existence 
of localization transition for eigenfunctions of the Schroedinger equation with disorder. They 
argued that ImS takes a finite value with non-zero probability. This means that in the 
absence of an external field the system 



Pk + «fcV^ 



i? 3 
i^k rx 



= 0, (51) 

R=Rfe R; 



has delocalized solutions.^ In other words, by analogy to the Schoedinger equation in the 
tight-binding approximation, the eigenstates {p^} of the system (|5T| ) extend throughout the 
entire volume. 

Fleishman and Andersonlll considered a three-dimensional system and neglected the 
Coulomb correlations. Note that if one rewrites the system (|5T|) in the momentum repre- 
sentation and averages it by factorizing the average of the product ap (mean-field) then one 
arrives at the plasmon equation ([13|) . As it was demonstrated in Section ||, this equation 
has no propagating solutions in the absence of Coulomb correlations. Our central point is 
that with the Coulomb correlations the delocalized solution exists even at the mean-field 
level. This solution, that is surface plasmon, is specific for the two-dimensional system and 
is characterized by the dispersion law u(q). 

An interesting question that could be addressed within the same approach is how the 
Coulomb correlations modify the a.c. Hall conductivity a xy (co) of the Anderson insulator.illl 
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FIGURES 

FIG. 1. (a) Diagrammatic presentation of the density-density correlation function H(u>, q); bub- 
ble (ij) stands for Pij(uj) while dashed line corresponds to the Coulomb interaction v(q). Each 
vertex contains matrix element (i|e* qr |j). (b) The block connecting two bubbles (ij) is replaced by 
a plasmon line, (c) The renormalization of the bubble (ij) caused by interaction with plasmons. 

FIG. 2. The real (a) and imaginary (b) parts of conductance is plotted as a function of dimen- 
sionless frequency hu>/2Io for U/2Iq = (solid line), U/2Iq = 3 (long-dashed line), and U/2Iq = 8 
(dashed line). 
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